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Abstract 



The extra particles provided by Supersymmetry(SUSY) appear as natural candidates for ex- 
otica such as the missing Dark Matter of the Universe. 

The particle candidates for Dark Matter, and the basic elements of the minimal supersym- 
metric extension of the Standard Model are reviewed. The needed elements to compute the 
neutralino counting rate are given, considering new Yukawa loop corrections to the neutralino- 
neutralino-Higgs boson couplings when the neutralino is a Bino-like state. 
Including the new loop corrections in the Bino limit, the counting rate is improved in a 20% 
for several values of the soft breaking supersymmetric parameters. [] 



1 Dissertation at the Diploma Course in High Energy Physics. ICTP. Trieste. Italy. December, 2000. 
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Chapter 1 
Introduction 



The Dark Matter (DM) problem is one of the most interesting problems in modern physics. 
What is the Dark Matter composition? What is the DM energy density present in our Uni- 
verse? What are the DM Baryonic and Non-Baryonic components? These are questions 
without a clear answer. 

From several estimates it is known that the DM constitutes the dominant mass density in our 
Universe, and hence is of vital importance for predicting the future of the Universe. 
The DM problem establishes a very nice connection between Cosmology and Particle Physics, 
given the possibility of predicting new particle candidates using the well known Standard 
Model and its extensions. In the Standard Model, the possible DM candidate is the neutrino, 
while in the Minimal Supersymmetric Extension of the Standard Model (MSSM), the lightest 
supersymmetric particle (neutralino or sneutrino) appears as possible candidate. Many ex- 
periments are looking for DM candidates 0. If one these experiments finds DM candidates, 
it will provide a very important hint to clarify some fundamental ideas, such as the presence 
of supersymmetry at the weak scale. 

As the fundamental theory to break Supersymmetry (SUSY) spontaneously is not known, the 
prediction of SUSY models such as the MSSM depends of many free parameters introduced to 
break SUSY. For example in the MSSM the soft breaking terms are introduced and therefore 
we take many free parameters entering in the interactions and the mass spectrums in the 
model. In the MSSM, the neutralino considered in many cases the ligthest supersymmetric 
particle is defined as the combination of the neutral gauginos. Different limits can be obtained 
for this particle, depending of the relation between the soft breaking parameters entering in 
its mass matrix. In our case we will consider the neutralino in the Bino-limit as the Dark 
Mattter particle candidate. 

In the present dissertation we will review the different particle candidates for Non-Baryonic 
Dark Matter: the Light Neutrinos, the Weakly Interacting Massive Particles (WIMPs) and 
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the Axions, and the basic elements of the Minimal Supersymmetric Extension of the Standard 
Model, such as the relation between particles and their superpartners, the general structure 
of the MSSM lagrangian including the soft breaking terms, the mass spectrums of Squarks, 
Neutralinos and Higgs bosons, and the R-symmetry, which is used to define the neutralino as 
a stable particle. 

In the Chapter 4 we introduce the different effective interactions to compute the neutralino 
counting rate when the neutralino is a Bino-like state. In this case with the objective of im- 
prove the neutralino counting rate we consider new Yukawa loop corrections to the neutralino- 
neutralino-Higgs boson couplings. Numerical results of the counting rate and the relic density 
are showed for several values of the soft supersymmetric parameters, where we can appreciate 
the increment of the counting rate when the new loop corrections are taking into account. 



Chapter 2 

Particle Candidates for Dark Matter 



2.1 Introduction. 

The phrase Dark Matter signifies matter whose existence has been inferred only through its 
gravitational effects. Evidence for DM comes from different cosmological scales, from galatic 
scales of several kiloparsecs and clusters of galaxies (Megaparsecs) to global scales of hundreds 
of Mpc's. 

The most well known evidence of DM comes from measurements of rotational velocities v TO t(r) 
of spiral galaxies. In this case at distances r greater than the extent of the visible mass, one 
expects that v TOt (r) oc r" 1 / 2 , however it remains constant. Therefore from this discrepancy, 
the existence of dark haloes is inferred to explain this effect. 

From different analyses in cosmology it is well known that the evolution of our Universe 
depends on the matter-energy present. Usually one uses the quantity Q = p/p c , where 
p is the average density of our Universe and p c the critical density. The critical density 
is given by p c = 1.9 x 1CT 29 h 2 g cm -3 , when the Hubble constant Ho is parametrised as 
h = #0/(100 km s _1 Mpc _1 ), the recent measurements of the Hubble constant give H = 
73 ± 6(stat) ± 8(syst) km s" 1 Mpc" 1 , so that 0.6 < h < 0.9. 

An estimate of the total matter density can be made from fluctuations in the Cosmic Mi- 
crowave Background Radiation (CMBR). The current data gives fixotai — 0.8 ± 0.2, in agree- 
ment with ^Totai — 1 predicted by models of cosmic inflation ||. This can be compared with 
the average density of visible matter in the Universe. 

The matter mass density may be derived from astronomical determinations of the average 
mass-to-light ratio M/L for various astrophysical objects (see for example ref. [BJ): 
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V M = ^ (2.1) 
from the M/L ~ (3 — 9)M Q /L Q ratio measured in stars, is possible estimate fi visible : 

0.002/T 1 < fi visib i < 0.006/T 1 (2.2) 

with upper value fusible — 0.012 

Also through measurements of X-ray emission in clusters of galaxies we can estimate the 
total matter density 0: 



'clusters 



0.4 ±0.1 (2.3) 



On the other hand, from analyses of Big Bang nucleosynthesis, comparing the observed abun- 
dances of H, He and Li, and taking into account the chemical evolution of the Universe due 
to stellar burning, one finds f2j: 

0.008 < fifaaryons/i 2 < 0.04 (2.4) 

where the upper limit in this case is ^baryons — 0.1. 

Note that from the previous evaluation we can conclude that the major part of matter in 
the Universe is not visible and some of this Dark Matter is not baryonic. Therefore we can 
define two categories: Baryonic Dark Matter, composed of baryons which are not seen, and 
Non-Baryonic Dark Matter, composed of massive neutrinos, or elementary particles which are 
as yet undiscovered, where the particles which comprise non-baryonic DM must have survived 
from the Big Bang, and therefore must be stable or have lifetimes in excess of the current age 
of the Universe. 

Also a very important point is that for purposes of galaxy formation models, the Non-Baryonic 
Dark Matter is classified as hot or cold, depending on whether the DM particles were moving 
at relativistic or non-relativistic speeds at the time galaxies could just start to form. 
Many candidates to describe the DM are known, ranging from the axions with m < 10 _2 eV = 
9 x 1O" 69 M , to black holes of mass M = 1O 4 M . 

The main baryonic candidates are massive compact halo objects (MACHOs), these include 
brown dwarfs, Jupiters, stellar black holes remnants, white dwarfs, and neutron stars. 
As for non-baryonic candidates, we have light neutrinos as a candidate to describe the hot 
component, and to describe the cold part, we have axions and weakly interacting massive 
particles (WIMPs). The axion is motivated as a possible solution to the strong CP problem, 
while WIMPs are stable particles, which arise in extensions of the Standard Model, such as 
heavy fourth generation Dirac and Majorana neutrinos, and in SUSY models the lightest 
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neutralino and sneutrino. 

In the next sections we will describe some of the particles candidates for Non-baryonic DM. 



2.2 Neutrinos. 

In units of critical density, the cosmic mass density for massive neutrinos is estimated as: 

n vh * = y 

j-i 93 eV 

The observed age of the Universe together with the measured expansion rate yields flDAih 2 < 
0.4 so that for any of the three families M: 

m v < 40 eV (2.5) 

This limit is probably the most important astrophysical contribution to neutrino physics. In 
this case, neutrinos with a mass 4 eV < m^ < 40 eV could represent all the non-baryonic 
Dark Matter. 

When we consider the low mass neutrinos as particle DM candidates, two problems are 
present: from the perspective of structure formation, neutrinos as hot dark matter would 
require topological defects such as cosmic strings to form small structures such as galaxies. 
The second problem is well known from the phase space constraint, for example if the neutrinos 
are the DM of dwarf galaxies, m u must be few hundreds of eV, which violates the upper limit 
from the overall cosmic mass density of about 40 eV 0]. Therefore these neutrinos can not 
describe the haloes of dwarf galaxies. However there exists the possibility to consider other 
candidates, for example the neutralino [[[]. 



2.3 WIMPs. 

The weakly interacting massive particles (WIMPs) are particles with masses between lOGeV 
and a few TeV, and with cross section approximately weak strength. Their relic density is 
approximately given by: 

„ , 9 0.1 pb c . . 

fl wlMP h 2 ~ - p 2.6 
{cr a v} 

where c is the speed of light, a a is the total annihilation cross section of a pair of WIMPs 
into SM particles, v is the relative velocity between the two WIMPs in their cms system, and 
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(. . .) denotes thermal averaging. 

Since WIMPs annihilate with very roughly weak interaction strength, it is natural to assume 
that their interaction with normal matter is also approximately of this strength. This raises 
the hope of detecting relic WIMPs directly, by observing their scattering off nuclei in a 
detector. 

Perhaps the most obvious WIMP candidate is a heavy neutrino, however a SU(2) doublet 
neutrino will have a too small relic density, if its mass exceeds a few GeV. Also, one has to 
require the neutrino to be stable. It is not obvious why a massive neutrino should not be 
allowed to decay. However in SUSY models with exact R-parity, the lightest supersymmetric 
particle (LSP) is absolutely stable. In this case we have the sneutrino and the neutralino 
as candidates. However the sneutrinos have a large annihilation cross section, their masses 
would have to exceed several hundred GeV for them to make a good DM candidate, this 
would be uncomfortably heavy for the lightest superparticle. 

The neutralino still can make a good cold DM candidate |5[]. It is well known that for very 
reasonable SUSY parameters the relic density is 0.1 < ^~o^ 2 — 0.3, which is the currently 
preferred range to the total DM density. 

In the next chapter we will analyze some properties of the neutralinos, such as masses and 
the different physical states. 



2.4 Axions 

Axions are hypothetical pseudoscalar bosons, introduced to solve the strong CP-problem of 

QCD using the Peccei-Quinn Mechanism || 

The axion mass arises from its mixing with the 7r°: 

m a ~ 0.6 e y 1Q7 ^ eV (2.7) 

la 

In units of the cosmic critical density one finds for the axionic mass density: 

n a h* ~ 0.23 x lff^«(_ i ^_)i-«^F(e < ) (2.8) 

where Of is related to the average initial value of the axion field and F(9i) with F(0) = 1 
encapsules anharmonic corrections to the axion potential. 

With Q a h 2 ~ 0.3 x 2 ±1 the mass of dark matter axions is found to be m a = 30 — 1000/xeV, 
astrophysics limit very important to take into account for future experiments ||. 



Chapter 3 

The Minimal SUSY Extension of the 
SM (MSSM). 

3.1 Introduction. 

The Glashow-Weinberg-Salam Model or Standard Model (SM), describes with a very good 
precision all electroweak processes. It is based on gauge invariance under the symmetry group: 

G SM = SU(3) C x SU(2) L x U(1) Y (3.1) 

and its partial spontaneous symmetry breaking. In table 1 we show all its constituents, the 
elementary fermions (quarks and leptons), the scalar Higgs boson, the gauge bosons, and their 
transformation properties under Gsm (using Q = T 3 + Y/2). The lagrangian of the Standard 
Model has the following form: 

^SM £ fermions ^gauge bosons ^scalars ^Yukawa (^'^) 

The explicit form of the SM lagrangian is well known, for our objectives we will write explicitly 
only the expression of the Yukawa interaction of the third generation: 

t-'Yukawa 

hit b) L <$>b R + ht(t b) L ia 2 <$>*t R + h T (P T t) l <$>t r + h.c (3.3) 

where hb, h t , h T , t, b, v T and r are the Yukawa couplings for each quark and lepton, the top 
quark, the bottom quark, the r neutrino and the r lepton, respectively (with R or L chiral- 
ity). Note that we can write these terms using only one scalar field <3>, which after electroweak 
symmetry breaking can generate mass for all the quarks and leptons in the SM. 
As was mentioned above, the standard model has an extremely economical Higgs sector, which 
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accounts for all the particle masses. Baryon (B) and Lepton (L) numbers are automatically 
conserved and it is an anomaly free Quantum Field Theory. However all is not perfect. At 
present there is no unambiguous direct experimental evidence for the Higgs boson. Moreover 
there are no answers to some fundamental questions such as: Could we explain the origin of 
parity and time reversal breaking? Could we unify the diferent gauge couplings? Could we 
explain the origin of the CP violation? 

Another fundamental problem in the Standard Model is in the Higgs sector. There is nothing 
to protect the mass of the Higgs boson from quadratic divergences. To solve this problem 
we can introduce a new symmetry called supersymmetry (SUSY). SUSY is a symmetry that 
transforms fermions into bosons and vice versa [71] . 

Motivations for introducing SUSY are the Haag-Lopuszanski-Sohnius Theorem (SUSY is the 
most general symmetry of the S-matrix), and cancellation of the tachyon states in String 
Theory. However the best (the more realistic) motivation is introduce SUSY at the weak 
scale to solve the problem related with the Higgs mass protection. 

Once SUSY is introduced we will have a superpartner for each particle in the Minimal Super- 
symmetry Standard Model, and one aditional Higgs and its superpartner are introduced [§]. 
Note that if SUSY is realized in Nature, it means that we have discovered only aproximately 
the fifty percent of all particles. 

From experiments (e.g at LEP) we know that the selectron, the superpartner of the electron 
must have a mass at least 90 GeV, more than 5 orden of magnitude above the electron mass. 
Therefore SUSY must be broken. This raises an additional fundamental problem: How can 
we break SUSY spontaneously? We know that when a symmetry is broken spontaneously we 
can keep the renormalizability and the unitarity of our model, but in this case we don't know 
the fundamental theory to break SUSY spontaneously. 

In the MSSM, SUSY is broken explicitly by introducing soft breaking terms in the lagrangian. 
These terms are called soft because the quadratic divergences (mass protection for the Higgs 
boson) are still canceled, but SUSY is broken ||. 

In the present chapter we will describe some properties of the MSSM, the mass spectrum of 
the Neutralinos, Squarks and Higgs bosons, and the main implications of the R-symmetry 
introduced in this model. 
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Quarks: 




4 

where i = 1, 2, 3 (colors) 
Leptons: 

en 
Scalars: 

-(?) 

Gauge bosons: 
G" wit/i a = 1..8 
wita 6 = 1..3 

Table 1. Standard Model Particles 



SU(3) C , SU(2) L ,U(1) Y 







(3r,2 L , 1/3) 


4 


t R 


(3c, U, 4/3) 


J. 
S R 


b R 


(3 c ,l L ,-2/3) 


V ) L 




(1c,2l,-1) 


Hr 




(lc,U,-2) 



(lc,2 L ,l) 

(8c,1l,0) 
(1c,3 l ,0) 
(1c,1l,0) 
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3.2 Particles and their Superpartners. 

In the MSSM, we add a superpartner for each SM particles in the same representation of 
the gauge group. In Table 2 we show the names of the superpartners. For example the 
superpartner of the top quark is called stop, of the photon the superpartner is the photino, 
and similarly for the other particles. However, in the Higgs sector we must add another new 
Higgs boson and its superpartner to write the Yukawa interactions needed to generate masses 
for all quarks and to obtain an anomaly free model. 



Table 2. Names of superpartners. 



Fermions -<=>- Sfermions 
(quarks, leptons) (squarks, sleptons) 

s = l/2 s = 

Gauge Bosons -<=>- Gauginos 
(W ± , Z, 7, gluons) (Wino, Zino, photino, gluinos) 

8=1 s = 1/2 

Higgs -<=>- Higgsinos 

s = s = l/2 

Usually we use the SUSY operators in the left-chiral representation, therefore it is con- 
venient to rewrite the SM particles (Table 1.) in the left-chiral representation and define 
the superpartners accordingly. This leads to the superfield formalism, which makes it easier 
to construct SUSY invariant lagrangians [[UJ . In this case we must introduce for each SM 
particle one superfield, which contains the SM particle, its superpartner and an auxiliary 
unphysical field. In Table 3 we show the third generation of the SM particles, their super- 
partners and the superfields needed to write our lagrangian. Note that we have an extended 
Higgs sector and the color index is omitted. 
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Table 3. Content of the MSSM. 



(only the third generation) 



Superfields 
Vector Superfields 



Bosonic Fields Fermionic Fields 



G 



SM 



G a 3 



Gn 



G l 



Chiral Superfields 
Leptons 



Quarks 

Q 

T 
B 
Higgs 

#1 
Ho 



G a n with a = 1..8 



W h u with b = 1..3 



B„ 



'R 



Ql 



t L 

b L 



'R 



J R 



m 

Hi 
HI 



G a 
W b 
B 



t L 

b L 



t c 



b?. 



Hi 
Hi 

m 



(8c,U,0) 
(1c,3l,0) 
(1c,1l,0) 



;ic,2 L ,-l) 
(Wl,2) 



(3c,2 l ,1/3) 

(3fc,l L ,-4/3) 
(3c, U, 2/3) 



'lc, 2l, -1) 
(lc,2 L ,l) 
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3.3 The General Structure of the MSSM Lagrangian. 

We can divide the lagrangian of the minimal extension of the Standard Model into two 
fundamental parts, the SUSY invariant and the Soft breaking term: 

Cmssm = C-susy + Csoft (3-4) 
In general we can write the SUSY invariant term as: 

£-SUSY = £-gauge + C-leptons + ^-quarks + ^-Higgs + J ^ 2 Q W + h.C (3.5) 

Having specified the content of the MSSM in Table 3 (including only the third generation of 
matter fields), we have to define the different terms of the lagrangian. The term C gaU ge has 
the following form: 

tgauge = \j d 2 9[2Tr(W 3 W 3 ) + 2Tr(W 2 W 2 ) + W^W 1 } (3.6) 

with: 

W 3 = -±DDex V (-G 3 ) D a exp(G 3 ) (3.7) 
Gs = E ^Gg (3.8) 

a=l Z 

W 2 = ~DDexp(-G 2 ) D a exp(G 2 ) (3.9) 
G2 = E^ (3-10) 

6=1 Z 

Wl = -^DDD a G 1 (3.11) 

where A a and o a are the Gell-Mann and Pauli matrices, respectively. The SUSY covariant 
derivatives D and D are used in the left chiral representation: 

D L = ^ + 2ia*0d ll (3.12) 

and 

BE=-| (3.13) 



and 
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The lagrangians for gauge interactions of leptons, quarks and the Higgs bosons are: 

beptons = J d 2 9 d 2 9 I) exp(2g 2 G 2 + 5l yGi)L+ 

+ J d 2 6 d 2 6 rt exp(2g 2 G 2 + gJjGi)r (3.14) 

C qU arks = J d 2 6d 2 6Q^ exp(2g 3 G 3 + 2g 2 G 2 + g 1 ^-G 1 )Q+ 



+ J d 2 9 d 2 9 Tt eMdiY Gl ~ 93^ a )*G a 3 )T+ 

+ J d 2 6 d 2 6 St expG^Gx - g 3 (X a )*G a 3 )B (3.15) 
= / d 2 6 d 2 6 Hi exp(2g 2 G 2 + g 1 ^G 1 )H 1 + 

+ J d 2 9 d 2 9 H\ exp(2g 2 G 2 + g 1 ^G 1 )H 2 (3.16) 

where g 3 ,g 2 ,gi are the SU(3), SU(2) and U(l) coupling constants, respectively. The Yj rep- 
resent the hypercharges of the different superfields. 

We can write the superpotential as the sum of two terms, W=Wr+Wnr- The first conserves 
lepton (L) and baryon (B) numbers: 

W R = eiji-iiHiHi + h T H[L j r + h b H{Q j B + h t H 3 2 Q i T] (3.17) 

where is the antisymmetric tensor, fx the Higgs mass parameter and h T , hb and h t are the 
different Yukawa couplings. The term Wnr, which explicitly breaks L and B numbers, is: 

W NR = eij [-iJ HiU + \VUt + XUQ j B] + \"BBT (3.18) 

In Wnr the first three terms break lepton number, while the last term breaks baryon number. 
If both B and L were broken, the proton would decay very rapidly, therefore the products 
A * A and A * A must be very small. In the last section of the present chapter we will analyze 
the R-symmetry related with the L and B numbers conservation and its implications. 
The soft breaking term is: 
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C s0 ft = mjl^l 2 + ml\H 2 \ 2 + m\ 2 {H x H 2 + H*H%)+ 



+Q L M~Ql + tnSmfctR + b R m\ b R + L L T M|L L + f^m 2 ~ R f^+ 



H 2 Q L (h t A t )t R - HxQ L {h b A h )b R - H 1 L L (h T A T )T R + 




(3.19) 



Note that in order to describe SUSY breaking we introduce many free parameters, and several 
terms have mass dimension less than 4 (super- renormalizable, but not SUSY invariant). Also 
the different mass terms considered lift the degeneracy between particles and their superpart- 



3.4 Neutralinos. 

Once SU(2) L x U(1) Y is broken in the MSSM, fields with different SU(2) L x U(1) Y quantum 
numbers can mix, if they have the same SU(3)c x U(l) em quantum numbers, and the same 
spin. 

The neutralinos are mixtures of the B, the neutral W and the two neutral Higgsinos. In gen- 
eral these states form four distinct Majorana fermions, which are eigenstates of the symmetric 
mass matrix [in the basis (B, W, Hf, H$)] ||: 



where sw = sin^ and Cw = cos 6^. 

Note that the masses of the neutralinos are determined by the values of four parameters M±, 

(H°) 

M 2 , jj and tan/3 = — = jtM- The neutralino mass term is given by: 



ners IITTI . 



M 



f M x 

M 2 

—Mz cos P sw Mz cos (3 cw 

\ Mz sin/3 sw —Mz sin/3 cw 



-Mz cos/3 sw Mz sin/3 sw \ 

Mz cos P cw —Mz sin P cw 

-n 

-li J 



(3.20) 
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- £™ ss = ~(V>°) T M </>° + h.c (3.21) 

If we define the physical states as Xi = Nijip®, the diagonal mass matrix is M D \".\/„.Y . 
In the limit Mi, M 2 , \i ^ Mz, we can diagonalize the mass matrix very easily pertubatively. 



In this case the eigenstates and their masses are |L2 
Bino-like: 

%i ^ B + N 13 H° + N U H% (3.22) 

where: 

M z sin6 l vy(Micos/? + y usin/3) 

N * ^zrwi (3 - 23) 



M z sin6 l H/(Misin/? + /xcos/?) 

Nu ~ q? + mi (3 - 24) 

and mass: 

m~o ~ Mi (3.25) 

Wino-like: 

X° 2 ^W + N 23 H° + N 2 JT$ (3.26) 

with 

_ M z cos 9 W (M 2 cos P + fxsmp) 

N23 ~ + Mf (3 ' 27) 



M z cos0iy(M 2 sin/3 + //cos/3) 
^ = J^Wi ( 3 - 28 ) 

and mass: 

m~o ~ M 2 (3.29) 

Higgsino-like: 

^^m-Hl NziB + NvW ,N u Hl 
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where 



and 



Si is: 



Higgsino-like: 



with 



M z sin0 w (sin/3 + cos/3) 
Nu ~ (~^ + M 1 ) (3 - 31) 



M z cos 6 W (sin (3 + cos (3) 
N " 2 ~ (J^Mi) (3 - 32) 



_ 5i cos26 w 
Nu ~ ~l + sin20 (3 - 33) 



m~o ~ /i(l + <Ji) (3.34) 



M|(l + sin2/3) sin 2 w cos 2 ^ 
5l " ^ '(^mJ+G^M^ (3 - 35) 



-L^ + + ^ ,3,6) 



Mzsin6W(-sin/3 + cos/3) 

^ " ( } 



= M z cos^(sin/3 + cos/3) 
42 ~ (/i + M 2 ) 1 • J 



iV 44 = (3.39) 
1 — sin 20 
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with 



m~o ~ -/i(l + 6 2 ) 



So is: 



M|(l - sin 2/3) sin 2 



w 



2/i 



+ 



cos 2 6 



w 



'(/i + Ml) (/i + M; 



(3.40) 



(3.41) 



Usually the unification of the gauge couplings is assumed, in this case the parameter M 2 ~ 
2M\. Therefore from this relation we can conclude that the Wino eigenstates never will 
correpond to the LSP particle. 

Now if Mz <C Mi < M2 < yU the LSP will be a nearly pure Bino eigenstate, while if 
Mz <C < Mi < M 2 the LSP is a linear combination of the Higgsino states. For more details 
of neutralino masses in the MSSM see ITBf . 



3.5 S quarks. 

After the electroweak symmetry breaking, several terms in the MSSM lagrangian contribute 
to the squarks mass matrices. Ignoring flavor mixing between sfermions, the mass matrices 
for the stop and sbottom squarks are ||: 
stop mass matrix: 



M4 



rn 



t + Ml + m\{\ - § sin 2 6 W ) cos 2/3 
Q 



-m t (A t + //cot/9) 



V 



sbottom mass matrix: 



-m t (A t + /icot (3) 
m 2 + m~ +| m| sin 2 6 W cos 2/3 



(3.42) 



M~ 



mj 



+ Ml- m|(| - \ sin 2 6 W ) cos 2(3 
Q 



-m^Ab + /ttan/9) 



— m^Af, + /xtan/J) m? + m2 

in the basis (/l,/r)- The physical states are defined as: 



\m 2 z sin 2 cos 2(3 



V 



(3.43) 
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fr = f L cos9 7 +f R sm9 r (3.44) 

and 

h = -f L sm9j+ f R cos6j (3.45) 
Note the contributions of the different soft breaking parameters in the mass matrices. 

3.6 The Higgs Bosons in the MSSM. 

As we have mentioned before, the existence of the Higgs bosons is the main motivation to 
introduce SUSY at the electroweak scale. In the MSSM the tree-level Higgs potential is given 
byi: 

Vm gg s = ™^|#i| 2 + m 2 H2 \H 2 \ 2 + m\ 2 (H x H 2 + h.c)+ (3.46) 
+ d±A(l Hl f _ |# 2 |2)2 + 9 l{HlH 2 f (3.47) 

o 2 
where m 2 H . = mf + |/i| 2 with i = 1, 2 

Note that in this equation the strength of the quartic interactions is determined by the gauge 
couplings. 

After electroweak symmetry breaking, three of the eight degrees of freedom contained in the 
two Higgs boson doublets get eaten by the W ± and Z gauge bosons. The five physical degrees 
of freedom that remain form a neutral pseudoscalar boson A , two neutral scalar Higgs bosons 
h° and H°, and two charged Higgs bosons H + and H~ . The physical pseudoscalar Higgs boson 
A° is a linear combination of the imaginary parts of H® and H®, which have the mass matrix [ 

m the basis (-^p, —j^U- 

M 2 = ( -m\ 2 tan/3 -m\ 2 \ 

1 I -"^12 ~ m l2 COt (3 J V ■ ) 

m\ = trM] = -2m\ 2 / sin(2/?) (3.49) 

The other neutral Higgs bosons are mixtures of the real parts of Hf and H 2 , with tree-level 
mass matrix [ ^ , ^ \: 

M 2 = ( -m\ 2 tan (3 + m| cos 2 (3 m\ 2 - \m\ sin 2/3 \ 

R \ m 12 ~ 2" m l Sin ^ - m 12 cot + m Z sin2 P J ' 
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In this case the eigenvalues are: 



m H°,h° = ot m A° +m z ± V ( m A0 +"i|) 2 - 4m^ m|cos 2 (2/3)] (3.51) 



From this equation we can see that at tree level, the MSSM predict that m^o < mz, unfortu- 
nately in this case when we consider one-loop corrections, the mass of the light Higgs boson 
is modified significantly. For example assuming that the stop masses do not exceed 1 TeV, 
nT>h° < 130 GeV 0| . In this example we can see the importance of loop corrections in the 



computation of physical quantities. 



3.7 R-symmetry and its Implications. 

In the Standard Model the conservation of Baryon (B) and Lepton (L) number is automatic, 
this is an accidental consequence of the gauge group and matter content. In the MSSM, as we 
showed in the second section of this chapter, we can separate the most general gauge invariant 
superpotential basically into two fundamental parts, where the first term conserves B and L, 
while the second breaks these symmetries. 

In the MSSM, B and L conservation can be related to a new discrete symmetry, which can be 
used to classify the two kinds of contributions to the superpotential. This symmetry is the 
R-parity, defined as: 

R = (_!)3(B-L) ( 3>52 ) 

Quark and lepton supermultiplets have R — — 1, while the Higgs and gauge supermultiplets 
have R = +1. The symmetry principle in this case will be that a term in the lagrangian is 
allowed only if the product of the R parities is equal 1. 

The conservation of R-parity as defined in equation fl3.52j ), together with spin conservation, 
also implies the conservation of another discrete symmetry, defined such that it will be +1 
for the SM particles and -1 for all the sparticles: 

R s = (_!)3(B-L)+2S (3>53) 

Now if we impose the conservation of Rs, we will have some important phenomeno logical 
consequences: 

• The lightest particle with Rs = — 1, called the lightest supersymmetric particle (LSP), 
must be stable. 

• Each sparticle other than the LSP must decay into a state with an odd number of LSPs. 
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Chapter 4 



Applications of SUSY Theories to the 
Description of DM. 

4.1 Introduction. 

Two methods have been proposed to search for relic neutralinos Jl|. The first is the direct 
detection, observing their scattering off nuclei in a detector. The second is to look for signals 
for ongoing neutralino annihilation. In the second case the idea is that the neutralino can lose 
energy in collisions with nuclei and can then be trapped by the gravitational field of celestial 
bodies, then eventually will become sufficiently concentrated near the center of these bodies 
to annihilate with significant rate. As is well known to analyze these possibilities, it is very 
important to know very well the neutralino-nucleon scattering cross section. 
In this chapter, knowing the basic elements of the SUSY extension of SM, such as interactions 
and spectrum, and assuming the neutralino as the lightest supersymmetric particle || , we will 
compute the counting rate in the limit where the neutralino is a Bino-like state. In this case 
we will include the one-loop corrections to the neutralino-neutralino-Higgs bosons couplings. 

4.2 Yukawa One-Loop Corrections to the Bino-Bino- 
Higgs Boson Couplings. 

As we mentioned in the previous chapter, we can obtain different limits for the neutralino 
physical states. One of the more interesting limits is the bino limit. In this case it is possible 
to compute with good precision the dark matter density contribution needed to describe the 
dark haloes in the spiral galaxies. 
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In this section we will analyze the LSP-LSP-Higgs bosons couplings in the Bino limit at one- 
loop level, with the objective of computing the counting rate in direct experiments. 
At tree level these couplings vanish in the pure Bino limit, therefore we will consider the 
one-loop correction to the elastic cross section. 



Using the Feynman Rules |Lq| , the tree-level couplings are given by: 



^% e o e ~o~ = z(sin a N 14 - cos a N 13 )(g 2 N 12 - 9l N n ) (4.1) 

Ko^o = *(sin a N 13 + cos a N u ){g 2 N 12 - 9l N u ) (4.2) 

^o-ojo = 7s(sin/5 N 13 - cos/3 N u )(g 2 N 12 - 9l N n ) (4.3) 

Note that for the pure Bino state (only Nu ^ 0) all these couplings are zero at tree level. 
Here the Higgs mixing angle a can be obtained from the following expressions: 

sin 2a = - sin 2/3 "f + ™f (4.4) 

and 







m 2 H o 


-m 2 h0 






m 2 A0 


— m 2 z o 



cos 2a = - cos 2(3^f f- (4.5) 

TTlljO m h° 

At one-loop level we will consider contributions which contain Yukawa interactions in the 
quark-squark loops. We therefore have to include heavy quarks and their superpartners inside 
the loops (see Fig |5.1| ). In this case the expressions for these couplings are given by: 
- Scalar Higgs bosons h° and if : 



, T 32 r(h°,H°) . , r(h°,H% A r 2 , , ( , , 

t6T x°x°(h°,H°) = ~^ ht6t + h ^l '\N n (4.6) 

where: 



i=l 



+2m~ m q J2(al + &|)[C (g~) - 2Cf (g)]} 



i=i 
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+ E G ~-i H "W^. - W~)C (q t , qj) + 2m~ (a~a~ + ftj^JCf («, $)] (4.7) 



-Pseudoscalar Hide's boson 

r 



^ F vOvO,o = + (4.8) 



with: 

Sf = Cf {£(a| - &|)(-m| + + mjL)tf (ft) + 2m~ m g £(a| + &|)C (ft)} 
«=i i=i 

2 

+ E C4°[m g (a~ 6~ ~ a ^) G o(Qi, Qj) ~ ^-oi^bqj + % b qJC±(qi, qj)} (4.9) 

Here c^ h ° ,H ° ,A °\ are the quark-quark-Higgs bosons couplings, while C-' 1 . ,H ,A ^ represents the 
q% — qj — Higgs bosons couplings: i, j labels the squark mass eigenstates. 
The interaction quark-squark-neutralino is written as: 

C ?m = ?K + ^)X% + h.c (4.10) 
We are using the following notation for the loop integrals |16|] (see appendix A): 

Co(gi) = Co(h, h,mq,m~,m q ) (4.11) 

CoiQuQj) = C (k 1 ,k 2 ,m~,m q ,m~) (4.12) 

Ct{qi) = Cf(-hi, k 2 , m q , rn q , m~) (4.13) 

CKq-hcfj) =C?(k 1 ,-k 2 ,rn~,m~,m q ) (4.14) 

Ci(.Qi,Qj) = c i(ki,-k 2 ,m~,m~,m q ) (4.15) 

where ki represent the neutralino momentum. As already noted only consider the contribution 
of heavy quarks and their superpartners in the loop (only the third generation quarks). 
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In the case of the top quark, using the Feynman rules given in [15] the parameters are 
Yukawa coupling: 

h t 



92m t 



and 



t 2 



\f2m w sin (5 
g 2 ttm6 w 1 2 . 

—j=- (-- 008%+ -8111%) 

g 2 time w 1 2 
= t= — - sm fc- H — cos 6M 

# 2 tan0 w 1 2 . 
Or = t= COS Ur sm Or) 

^ 2 tan ^ 1 2 

Or = t= -Sin&r COS 6M 

* 2 V 6 * 3 * ; 

the couplings for the heavy Higgs boson H° are given by: 
t-t-H° 

sin a 



t-t-H c 



C 



V2 



C 



til 



^ m 2 w sin/?cos(a + 0) I 2 2 . 2 
-V2[ — ^ (-cos %--sm ^cos2% 



(4.16) 

(4.17) 
(4.18) 
(4.19) 
(4.20) 

(4.21) 



+ m t sin a sin 26^{/i cos a + A t sin a)] 

2 



(4.22) 



t22 



m 2 ^ sin j3 cos(a + /3) . 1 



m 4 cos 2 # 



w 



- sin 6*r H — sin 6 W cos 2%) 
2 1 3 



+ m t sin a + - sin 2#r(yU cos a + A t sin a)] 



(4.23) 



tl2 ~~ <21 



m 2 ^ sin /3 cos(a + (3) sin 2% 



4m* cos 2 6 



w 



— - cos 26*r(yU cos a + A t sin a)] 



(4.24) 
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In the case of the light Higgs boson the couplings are: 
t-t-hP 



t-t-h° 



Cf = (4.25) 



o = _^ m^sin/3sm( a + /3) 1 2 2 ^ 

tn L m t cos 2 w v 2 * 3 * y 



+ m t cos a + - sin sin a — A t cos a)] (4.26) 



/TTr ^ W /sin / 9sin(a + /?) 1 . 2 2 . 2 



+ m t cos a — ^ sin 26^(/i sin a — A t cos a)] (4.27) 



C h° = C h° = ^/2\ m w sin ft sin ( a + ft) sin 2 % 
* 12 * 21 ' " 4m t cos 2 W 



+ - cos 2^{// sin a — A t cosai)] (4.28) 



1 

2 

For the contribution of the pseudoscalar Higgs boson A , the couplings are: 
t-t-A° 



yA o .cos (3 



t-t-A° 



Cf = -i^r- (4.29) 



C £ = - C 7n = " A * cot ^ ( 43 °) 

C& = C& = (4.31) 



We now turn to the sbottom coupling. 
The Yukawa coupling is: 
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h = »Hi (4.32) 



\/2mw cos (3 

and 



a,~ = = — — cosfc'r sm9r) 4.33 

as= 9^ ( l sin9 ._l cos( y (4 34) 

# 2 tan0 w 1 1 . 

Or = f= — ( — cos Or + - sin Or) 4.35 

g 2 tanO w 1 . 1 

b ~ 2 = 72 ( 6 Sm ^ + 3 C ° S ^ ( ) 



The 6 — 6 — H° coupling is given by: 

V2 



Cf = (4.37) 



and the couplings b — b — H° are: 

^ w o /TTr^w cos/9cos(o; + /3) . 1 9 ^ 1 . <>„ 

= _V2 _^ — cos 2 % + - sin 2 ^ cos 2% 

611 1 m b cos 2 fl^ v 2 b 3 b; 

+ m;, cos a — - sin 2^-(/i sin a + Aj, cos a)] (4.38) 
r H ° - Kr ™ 2 w cos f3 cosja + f3) 1 2 1 2 pn , 9 ^ 



+ m& cos cH — sin 20r(/j, sin a + cos a)] (4.39) 
2 



_ _ ^r^w cos (3 cos(a + (3) sin 2^ 

^19 ~~ ^91 ~~ V 



Am b cos 2 



— i cos 26^(/x sin a + A b cos a)] (4.40) 
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In the case of the light Higgs boson h° the couplings are: 
b-b-h° 



b-b-h° 



Cf° = (4.41) 

o = ^ m^cos/3sin(a + /3) 1 1 
611 L m b cos 2 6 w v 2 6 3 6; 

+ m b sin a + ^ sin 29^(fi cos a — A b sin a)] (4.42) 



Wi° /- r mL cos/3sin(o; + /3) , 1 . 2 1 . 2 . 

+ rrib sin a sin 26^(/i cos a — A b sin a)] (4.43) 



b21 [ 4m b cos 2 9 w 



+ i cos 2^(yU cos a — A fe sina)] (4.44) 



The couplings of the pseudoscalar Higgs boson A : 
b-b-A° 



Cf = ^ (4.45) 



b - b - A° 

CA° = _ ^° = ^ 

^12 ^612 



Qf^-^^^^ + Atan/?) (4.46) 
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4.3 Effective Interactions. 

The neutralino-nucleus elastic cross section is of fundamental importance. It determines the 
detection rate in direct and indirect detection experiments. 

In general the WIMP-nucleus elastic cross section depends fundamentally on the WIMP- 
quark interaction strength, on the distribution of quarks in the nucleon and the distribution 
of nucleons in the nucleus 0. 

To compute the WIMP-nuclei interaction, we must first compute the interactions of WIMPs 
with quarks and gluons, after using the matrix elements of the quark and gluons operators in 
a nucleon state we compute the interaction with nucleons and finally using the nuclear wave 
functions, it is possible to compute the matrix element for the WIMP-nucleus cross section. 
In the non-relativistic limit we have only two types of neutralino-quark interaction, the spin- 
spin interaction and the scalar interaction. In the case of spin-spin interaction, the neutralino 
couples to the spin of the nucleus, and in the case of the scalar interaction the neutralino 
couples to the mass of the nucleus. 

To start with the analysis of the neutralino-nucleus interaction, we can discuss the effective 
lagrangian describing the neutralino-quark interaction. In this case we have three classes of 
contributions, the exchange of a Z° or Higgs boson and squark exchange (Fig [5.2| ). 



From the Z° and squark exchange we have spin-dependent contributions, while the Higgs and 
squark exchange contribute to the spin-independent interactions. 
The spin-dependent contribution is |T^| : 



L% = xV75X°g7„(c 9 + d ql ,)q (4.48) 



where: 



-~ £ 2 T k \ * + SrO R iT, q - 2e q sin 2 6 W ) (4.49) 

1 2 a 2 ~ + b 2 ~ al D 

d q = 7 E —2 f so - l^rO R T 3q (4.50) 

The second term in c q and d q represents the Z° exchange, g% is the SU(2) gauge coupling, 9w 
the weak mixing angle, T 3q = ±1/2 and e q the weak isospin and electric charge of the quark 
q, and R = ^(A^ 2 3 — iV 2 4 ) describes the Z°XiXi coupling. 

To determine the spin-dependent neutralino-nucleon interaction, we must compute the matrix 
element: 

(N\q lfll5 q\N) = 2s„AgW (4.51) 
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Here s M is the spin vector of the nucleon and Aq^ denotes the second moment of the polarized 
quark density. The Aq^ can be extracted from analyses of polarized deep-inelastic lepton- 
nucleon scattering Here will use the values: 

Ad^ = A«W = 0.77 (4.52) 
Au (n) = Ad (p) = -0.49 (4.53) 

As {n) = As {p) = -0.15 (4.54) 

Now knowing these values, we can write the effective spin-dependent interaction neutralino- 
nucleon as: _ 

C%{ n _ dep = x°Yl^NS^ N £ 2dM N) (4-55) 

q=u,d,s 

where $at denotes the nucleon wave function. Note that this interaction is considered at tree 
level. 

The spin-independent or scalar neutralino-quark interaction is usually the most important 
contribution. It is given by [17|: 



£ e /J n -in dep = f q X°X°qq + g^O^x^x ~ ^m^qqx^X ] (4.56) 



The leading twist-2 quark operator is: 



i 1 
°£l = 2 &7/Ag + qivd^q - -qd^qg^] (4.57) 



where 



1 2 a 



2 ~-b 2 ~ 2 C%CP 



and 



f q = - L Y—o f -T + m.V *\ q (4.58) 

" 4 ^ ml - (m~ + m q f q ^ m\, 1 } 



1 2 a 2 ~ + b 2 ~ 

S, = ~En — r ~ — s™ (4-59) 

4feK~-K~o + <> 22 



The couplings C~ and represent the couplings of the jth scalar Higgs boson to the LSP 
and to the quark q respectively. Including the one-loop corrections obtained in the previous 
section, we can write C~ as: 
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n {h ,H ,A ) _ r tree , e r (a an \ 

~~ (h°,H°,A a )x X + (/i°,H ,A )x x° v DU ' 

The different contributions to the LSP-gluons interaction are shown in (Fig |5.3|) . The effective 
interaction from these graphs is given by ||18|| : 

Cgiuons = b^x°G%G a ^ - (B 1D + B ls )7d,d u x°G^+ 

+ B^iid^ + id u ^)x°G (2 ^ (4.61) 

b = T~+T g + B D + B s - ~^B 2S - -*-(B 1D + B 1S ) (4.62) 



where: 



with 



b d = - b l) m M m g, ™v m ~o) (4-63) 

q,i 

a 2 1 

B s = -^--m~ Y(aZ + b 2 ~)I 2 (m~ m q , m~ ) (4.64) 
167r 2 8 x qi qi q x 

q,i 

g 2 i 

Bid = j^- XX a ! ~ h \) m ih( m v m n m ~o) ( 4 - 65 ) 

B 1S = T^o m ~0 + & |) J 4( m ^ "V m ~o) ( 4 - 66 ) 

^ = J— £(a| + ^K* m q , m~ ) (4.67) 
„ 2 ! 2 C<2> G i?' 



167r 12 t-[ m^ 

,2 , 2 C W 



lb\Y,z£- E C« (4.69) 



167T 3 r~f mtro , 

J=l ifV q=c,b.t 

Here T~ and T 9 are the Higgs contributions from squark and quark loops respectively, while all 
other contributions come from box and triangle diagrams involving quarks and squarks. All 
coefficients with subscript D are proportional to the difference a~ —b~, while the subscript S 

Qi Qi 
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i 2 ~+b 2 ~. 

We have used the leading twist-2 gluonic operator: 



indicates a contribution proportional to a~ + b~ 



The loop integrals as function of a Feynman parameter are given by: 



,1 x 2 -2x + 2/3 
h(m~,m q ,m x )= j dx — (4.71) 



l^i x(x 2 -2x + 2/3) fA „. 

I 2 (m~,m q ,m x ) = ^ dx^ — '—L (4.72) 



1 , x 2 (l — x 



2 



I 3 (m~,m q ,m x ) = dx — (4.73) 



1 3/ (l 

h(m~m q ,m x )= £ — (4.74) 



/ 1 , x(l — x)(2 — x) 
I 5 (m~,m q ,m x ) = [ dx-± ^ '- (4.75) 



with 



D = x 2 m 2 ~ + x{m~ — m 2 — m2 ~ ) + m l (4.76) 

Note that this effective lagrangian is speficied at a high energy scale, for example /xq — Tn~. 
With the objective of writing the effective neutralino-nucleon interaction we must evaluate 
the matrix elements of the quark and gluon operators in a nucleon state [l|]. 
For light quarks we write (q = u, d, s): 

(N\m q qq\N) = m N ff q ] (4.77) 

where for the neutron we have:/y„ = 0.023, ftf 1 ] — 0.034, f!pj = 0.14, and for the proton: 

= 0.019, f^ d = 0.041 and = 0.14. 
The heavy quark matrix elements can be computed perturbatively: 

(N\m Q QQ\N) = Jwi - £ f$\ = ^m N f TO (4.78) 

q=u,d,s 

We must also compute the matrix elements of the ( - 2 ^ fJ " u and G^ 2 ^ v operators, at zero mo- 
mentum transfer M: 
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m N (N\0^\N) = (pV - -m 2 N <T) [' dx x [q(x, /i 2 ) + q(x, /i 2 )] (4.79) 



m N (N\G^\N) = ( P y - -m 2 N <T) f dx x [G(x, /i 2 ) + G(x, /i 2 )] (4.80) 

4 jo 

where is the momentum of the nucleoli. The functions q(x, /i 2 ) and G(x, /i 2 ) are the quark 
and gluon densities in the nucleon at the scale [i. 

Using all the matrix element mentioned before, we can write the effective lagrangian for the 
scalar neutralino-nucleon interaction: 

*£L = f {N) 7 X°®N®N (4.81) 

where: 

f(N) fW m~ m q 2 f { q H) 



m N , m a 2 27 m a 

JV q=u,d,s 1 q=c,b,t 1 



- -f- £ 9* itf) ~ + ~f- G ^ 2 ) 93^) [B 2S + ^rn~ (B 1D + B lS )) (4.82) 

q=u,d,s,c,b " ~~ 

Here f^ H ' is the Higgs boson exchange contribution to the coefficient f q : 

, 2 g 2 (N 3 isina + cosa)(N 2 i — Nntan9w)h H o 
/f ) = £ —J * (4.83) 



i=l 4m ^o 



where h H o qq : 



h>H°tt = 


m t sin a 
sin (3 


(4.84) 


hh°tt = ■ 


m t cos a 
raw sin /? 


(4.85) 


hH°bb = 


cos a 
mw cos /3 


(4.86) 


hh%b = 


mi, sin a 
mw cos /3 


(4.87) 



The new loop-corrections computed appear in the Higgs exchange contribution to f q (eq |4.58| ), 
as well as in f( H > and the Higgs exchange contributions T~to b (equation |4.62|) . 
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4.4 The Counting Rate. 

The most direct possibility to detect DM particles is look for their scattering off the nuclei 
in a low-background detector. In general the counting rate is defined as events per days per 
kilogram of target, depending on the density and velocity distribution of dark matter particles 
in the galactic halo near the Earth. In their velocity distribution, one should take into account 
the motion of the Sun and Earth, which increase the total rate and give a yearly modulation 
to the differential rate. It will be important as method of distinguishing signal and noise if 



many events are found. More details about nuclear physics of DM can be found in |19] and 
i ' 

In general the differential detection rate can be written as [|i~7f|: 



N T ^ fdvf(v)4^(v,E R ) (4.f 



dE R m~ J dE_ 

where Nt = 6 - 02 * 1q26 is the number of the target nuclei per kg, p~ is the local neutralino 
matter density, v and f(v) denote the velocity and the velocity distribution function in the 
Earth frame, Er = mrV ^ € ° sg - is the nuclear recoil energy, where 9* is the scattering angle in 
the neutralino-nucleus center of mass frame, is the nuclear mass and m r the neutralino- 
nucleus reduced mass. 

Considering a Maxwell distribution for the galatic neutralinos, a Gaussian nuclear form factor 
and integrating over all the possible incoming velocities and deposited energies, the counting 
rate is: 

R - °Z l-SxlO^e^ P ? 

m ~o m N kg days 0.3 GeV/ cm 3 320 km/ s 

where is the average neutralino velocity, which results from the integration over the 
Maxwellian velocity distribution, a is the elastic neutralino-nucleus cross section at zero 
momentum transfer, it is given by: 

4m 2 ~m% 

a = -{[Zf& + (A — Z)/W] a + 4A 2 J(J + d q Aq( N) ?} (4-90) 

7T { m ~Q + m N) uAs 

We will consider as target the isotope 73 of Ge and A 2 J( J + 1) = 0.065 for this nucleus, where 
J is the total spin of the nucleus and A the nucleonic matrix element. 



0.573 r exp[-B/(l + g)]er/[l/^(l + g)] 

vtt^ } (4 - 91) 
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describes the suppression due to the Gaussian nuclear form factor considered. erf(. . .) is the 
error function and 



m 2 ~m 2 N o 

B = £ ^r 2 v~ (4.92) 

with V~ = f^. 

For the spin-independent interaction, we will use r spin _ indep = (0.3 + 0.89A 1 / 3 ) fm, while for 
spin- dependent interaction r spin _dc P = 1-25 r spin _ indep . 

Before showing some numerical results we must mention different aspects taken into account: 

• The neutralino is considered the lightest supersymmetric particle. It means that we 
must compute all the sparticle masses and impose this condition. 

• We will take into account the experimental constraints from SUSY search experiments, 
given in reference @: m~ > 90 GeV, m~ > 89 GeV and m~ > 43 GeV. The experi- 

Xi 

mental contraints on the Higgs search also play an important role. 

• Also we will show the values for the neutralino density, to estimate the range where Qh 2 
has an interesting value for cosmology. 

The relic density is computed, taking into account the annihilation of neutralinos in: 
two scalar neutral Higgs bosons, one scalar plus one pseudoscalar Higgs boson, / plus /, 
with / being an Standard Model fermion, Z plus Higgs bosons, W~ W + , Z Z, W + H~ , 
H~H + and the co- annihilation with charginos and sfermions. 

In our case, we are interested in the Bino limit, where M\ < M2, fi- In the soft breaking 
term, we will consider A t = Ab = A and Mq = m t ~ = m fc ~ for simplicity. 
Knowning the details related with the counting rate and the experimental contraints for 
SUSY particles, we will show some numerical examples, where we compute the counting rate 
including and excluding the new loop corrections to the LSP-LSP-Higgs boson coupling. 
For these calculations, we will use: = 400 GeV, m^o = 500 GeV, the Bino mass Mj = 
200 GeV, the Wino mass M 2 = 400 GeV and tan/5 = 5. The parameter \i is varied from 
±500 GeV to ±1050 GeV in general. For A we consider the values ±2.5, ±2.6 and ±2.7 

M Q- 

Note that in this case, the mass spectrum for supersymmetric particles is above the SM 
particles mass spectrum. It means that in this case we have not overlapping of these mass 
spectrums and all sparticle mass bounds are easily satisfied. 
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The first example is shown in Fig |5.4| and Fig |5.5| . In this case we show the values for the rate 
and the relic density respectively, when A = 2.5 and —1050 < fi < —500. The counting rate 
is incremented when the loop-corrections are considered, where the bigger difference is of a 
20 % in the interval —1050 < fi < —900. Analyzing the relic density, we have a significant 
value in the interval —1050 < \i < —950, and when the loop-corrections are included, it 
does not change appreciably. Note that the mass of the lighter stop squark increases with 
increasing \x. The upper bound on /i is set by the requirement that the neutralino must be the 
lightest SUSY particle. For the upper bound, the relic density is greatly reduced by x? _ ti 
co-annihilation . 

For the opposite choise of signs A = —2.5 and 500 < \i < 1050 (see Fig [5]6] and Fig |5.7|) , the 
counting rate does not change. The relic density for these values has a significant value in the 
interval 900 < fx < 1050 and it does not change, when the loop corrections are considered. 
Note that now m t ~ increases with increasing /i and the effect of ti — xl co-annihilation is again 
evident. 

In Fig and Fig |5.9j , the counting rate and the relic density for A = 2.6 and —1050 < 
fx < —700 are shown. The counting rate is incremented by a 16-20 % if the loop-corrections 
are included, while the relic density is quite small in this interval, equal to Qh 2 ~ 0.1 for 
-1050 < fi < -1000. 

Considering the loop-corrections, in the calculation of R and Qh 2 for A = —2.6 and 700 < 
/i < 1050 (Fig |5.10| and Fig [5.11| ), we see that R is increased and flh 2 remain equal, but it is 
very small in this case. 

In the last example, we consider \A\ = 2.7 and 900 < < 1050, taking into account two 
different cases for the different signs of these parameters. When A = 2.7 and fi is varied from 
-1050 to -900 (Fig [5. 12| and Fig |5.13| ), the inclusion of loop corrections increases R, while the 
relic density is very small in this interval. Changing the sign of these parameters(Figs |5.14 
and |5.15|) , R is increased and Qh 2 remains equal. 

We see that for the different values and signs considered of A and \i the Rate is increased, 
but the size of the correction is modest (~ 20%). 

The size of the corrections increases with decreasing m^. However, one has to respect the 
constraint m t ~ > ^~o- Moreover, for ~ m~ the predicted LSP relic density is small, due 

to co- annihilation. It should be noted, that in our calculation only include thermal relics, 
which were in thermal equilibrium after inflation. There might be significant contribution to 
Qh 2 , from non-thermal relics [^T[ 

In this way, we conclude that for the predictions of new experiments to search for SUSY Dark 
Matter candidates, might eventually be necessary to include loop corrections to the counting 
rate. However, these corrections are only important when the gaugino-higgsino mixing is 
small, i.e when the LSP is nearly pure Bino. In this case the predicted counting rate is 
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unfortunately several orders of magnitude below the sensitivity of present experiments. 



Chapter 5 
Conclusions. 



The different Non-Baryonic Dark Matter particle candidates and the basic elements of the 
Minimal Supersymmetric Extension of the Standard Model needed for our purposes were re- 
viewed. 

To improve the calculation of the neutralino counting rate were introduced new Yukawa loop 
corrections to the coupling neutralino- neutralino-Higgs boson in the Bino-limit. 
We have shown several examples to analyze the contributions of these corrections on the 
counting rate, considering the neutralino mass Mi = 200 GeV. The difference in the neu- 
tralino counting rate when the loop corrections are considered can reach 20%. Since these 
corrections are significant only when the counting rate is small, it is not important for present 
or near-future experiments. However, it plays an important role when one is trying to deter- 
mine the minimal realistic counting rate in the MSSM |f?2]j . 

The interesting interval for the relic density was taken into account to define when the loop 
corrections are important, the supersymmetric parameters were chosen according to the ex- 
periment constraints given at LEP2 ||. 

In general we conclude that it will be interesting to consider these Yukawa one-loop corrections 
in the different limits of neutralino physical states, and study the intervals of the different 
supersymmetric parameters which would provide the possibility of improving the counting 
rate appreciatly. 
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CHAPTER 5. CONCLUSIONS. 
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APPENDIX A. 



The loop integrals fll6 



with 



C Q (p 1 ,p 2 ,M 1 ,M 2 ,M 3 ) = (2tt//) 



4— n 



d n l 1 
i^DxD 2 D 3 



D l = l 2 - Ml + ie 



D 2 = {l + pif-Mj + ie 



Ds = (I+P1+P2) 2 ~ M$ + ie 



— , d n l 1 

C> 1 ,p 2 ,M 1 ,M 2 ,M3) = (2^) 4 - ,, ' 



in 2 D A D b D 6 



where: 



D 4 = (l-pi) 2 -Mf + ie 



D 5 = (l- p 2 y - Mi + ie 



D 6 = l 2 - Ml + ie 



C„ = C+(pi + p 2 )„ + G\ (pi - p 2 ). 



FIGURES. 




Figure 5.2: This figure shows the different contributions at tree level, to the neutralino-quark 
interaction. 




Figure 5.3: This figure shows the different contributions to the LSP-gluons interaction. 
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Figure 5.4: 




Figure 5.5: 
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Figure 5.7: 
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Figure 5.8: 




Figure 5.9: 



A= - 2.6 





Figure 5.11: 



A= 2.7 



0.000114 
0.000112 
0.00011 
0.000108 
0.000106 
0.000104 
0.000102 
1e-04 
9.8e-05 
9.6e-05 



9.4e-05 









loop-corrections 










no ioo^ 


j-correc 


lions -- 






































































































































-1060 -1040 -1020 -1000 



-980 



-960 -940 -920 -900 



Figure 5.12: 
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Figure 5.13: 



